I  : 

■7 

'  AD-A095  641  AIR 

'  LAM 

I  MAY 

1  UNCLASSIFIED  AFW 

FORCE  WEAPONS  LAB  KIRfLAND  AFB  NM 

MULTIBURST  CALCULATIONS  FOR  VARIOUS  ATTACK  S 
80  J  N  AUBREY#  H  J  ABEYTA#  W  E  GIFFORD 
.-TR-79-200 

CENARIOS. 

F/G  18/3 

U) 

NL 

ll 

1 

- 

l 

_ 

, 

Hi 

■■ 

■r 

■H 

■■ 

■i 

HI 

■■1 

■1 P*8* " 

f. 

■ _ 

AFWL-TR-79 


200^* 


f '  i 


fc?~»  r 

M‘v*  ,vN  «  &  **; 


LAMB  MULTIBURST  CALCULATIONS  FOR 
VARIOUS  ATTACK  SCENARIOS 


H.  J.  Abeyto 
Capt  W.  E.  Gifford,  III 
J.  W.  Aubrey 


May  1980 


Final  Report 


K.< 


AFWL-TR- 

79-200 


Approved  for  public  release;  distribution  unlimited. 


THIS  RESEARCH  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR 
AGENCY  UNDER  SUBTASK  Y99QAXSG677,  WORK  UNIT  BZO,  WORK 
UNIT  TITLE:  FREE-FIELD  AIRBLAST  ENVIRONMENT  DEFINITION. 


i 

'  Cxi 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305 

AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


81 


i  1 


AFWL-TR-79-200 


This  final  report  was  prepared  by  the  Air  Force  Weapons  Laboratory,  Kirtland 
Air  Force  Base,  New  Mexico,  under  Job  Order  WDNB4405.  Mr  James  W.  Aubrey  (NTESB) 
was  the  Laboratory  Project  Officer-in-Charge. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  a  definitely  related  Government  procurement  operation,  the 
Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever,  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded 
by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  persmission  to  manufacture,  use, 
or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  authored  by  an  employee  of  the  United  States  Government. 
The  United  States  Government  retains  a  nonexclusive,  royalty-free  license  to 
publish  or  reproduce  the  material  contained  herein,  or  allow  others  to  do  so, 
for  the  United  States  Government  purposes. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  avail¬ 
able  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


GARY  P.  GANONG 
Lt  Colonel,  USAF 

Chief,  Technology  and  Applications  Br 


FOR  THE  DIRECTOR 


Lt  Colonel,  USAF 

Chief,  Civil  Engineering  Research  Div 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  PAGE  (Whm  Data  Enl.r.dJ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I  if  _ 

•v.  (7g2 

■  lb 


LAMB  JJULTIBURST  Ct 
ATTACK  SCENARIOS* 


"U  '  nvn(  a/ 

0  W.^ubrey 
H.  J./\beyta 


ALCULATIONS  FOR  VARIOUS; 


RO.  H&’OfW  NUMBER 


W.  E^Gifford,  Ulj 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Air  Force  Weapons  Laboratory  (NTESB) 
Kirtland  Air  Force  Base,  NM  87117 


8.  CONTRACT  OR  GRANT  NUMBER/aJ 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


WDNB4405 

62704H 


Q&ZEs 


H.  CONTROLLING  OFFICE  NAME  AND  ADORESS  ■■  REPORT  QATpi... 

Director  fj/J  Mau  I 

Defense  Nuclear  Agency  v-^=wiFwui!lftEi^F t4ges 

Washington  D  C  20305 _ _ _ 1_36 _ _ _ 

1 4.  MONITORING  AGENCY  NAME  ft  ADDRESS^/  different  from  Controlling  Office)  15.  SECURITY  CLASS,  (of  thia  report) 


Air  Force  Weapons  Laboratory  (NTESB) 
Kirtland  Air  Force  Base,  NM  87117 

16.  DISTRIBUTION  STATEMENT  (a I  this  Repot!) 


UNCLASSIFIED 


1 5a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


Approved  for 


areer~dis  tTtbutTwrtml  fmtt 

y?9qA*P 


17.  DISTRIBUTION  STATEMENT  (0 


‘20,  If  HHerent  from  Report) 


L _ - J 


18.  SUPPLEMENTARY  NOTES 

This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask 
Y99QAXSG677,  Work  Unit  BZO,  Work  Unit  Title:  Free-Field  Airblast  Environment 
Definition. 

19.  KEY  WORDS  (Continue  on  reveree  aide  if  neceaamry  end  identify  by  block  number) 


LAMB  system 
Multiburst 
Overpressure 
Overpressure  Impulse 


Reentry  Vehicles 

Shock 

Nuclear 


20.  ABSTRACT  (Continue  on  reveree  aide  If  neceaamry  end  Identify  by  block  number) 

-^This  report  contains  the  results  of  33  Low  Altitude  Multiburst  (LAMB)  calcula¬ 
tions.  These  calculations  provide  overpressure  and  overpressure  impulse  time 
histories,  as  well  as  shock  front  velocities  for  12  multiburst  attack  scen¬ 
arios.  Yields  of  1 ,  3,  and  5  Mt  were  used  in  these  calculations,  with  shelter 
spacing  of  1000,  1500,  or  2000  m. 

/\ 


DO  1^*73  1473 


UNCLASSIFIED 

S E CURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Date  Entei 


AFWL-TR-79-200 


CONTENTS 

INTRODUCTION 
APPROACH 
RESULTS 
APPENDIX  A 
APPENDIX  B 


1/2 


AFWL-TR-79-200 

INTRODUCTION 

This  report  consists  of  33  LAMB  calculations  (Table  1)  and  a  description  of 
the  data.  It  was  written  to  provide  overpressure  and  overpressure  impl use-time 
histories  as  well  as  shock  front  velocities  for  various  multiburst  scenarios. 

The  Air  Force  Weapons  Laboratory  LAMB  code  was  used  because  of  its  fast  computa¬ 
tion  running  time  and  economy. 


TABLE  1 .  MULTIBURST  PROBLEM  MATRIX 


PEAKS 


Probl em 


RV/Shelter 

Attack 


Spacing 


1 

33.1012 

l/2(a) 

3  Mt 

0  -n 

1500  m 

3.5 

2.4 

2 

33.1122 

l/2(b) 

3  Mt 

0  IT; 

1500  m 

7.6 

1.9 

3 

33.1232 

l/2(c) 

3  Mt 

0  m 

1500  m 

3.6 

2.5 

4 

33.1142 

1  /2(d) 

3  Mt 

0  m 

1500  m 

7.6 

1.8 

5 

33.8001 

1/3 

3  Mt 

0  m 

1000  m 

18. 

2.8 

6 

33.8002 

1/3 

3  Mt 

0  m 

1500  m 

4.8 

1.3 

7 

33.8003 

1/3 

3  Mt 

0  m 

2000  m 

2.0 

0.7 

8 

53.8001 

1/3 

5  Mt 

0  m 

1000  m 

30. 

4.7 

9 

53.8002 

1/3 

5  Mt 

0  m 

1500  m 

7.8 

2.1 

10 

53.8003 

1/3 

5  Mt 

0  m 

2000  m 

3.2 

1.2 

11 

33.8012 

l/3(a) 

3  Mt 

0  m 

1500  m 

4.8 

1.4 

12 

53.8012 

1/3(a) 

5  Mt 

0  m 

1500  m 

7.8 

2.2 

13 

33.8032 

l/3(c) 

3  Mt 

0  m 

1500  m 

2.3 

1.4 

14 

53.8032 

1 /3(c) 

5  Mt 

0  m 

1500  m 

3.8 

2.5 

15 

33.1201 

1 /3(e) 

3  Mt 

0  m 

1000  m 

9.3 

4.6 

16 

33.1202 

1 /3(e) 

3  Mt 

0  m 

1500  m 

3.5 

2.1 

17 

33.1203 

l/3(e) 

3  Mt 

0  m 

2000  m 

1.7 

1.1 

18 

53.1201 

l/3(e) 

5  Mt 

0  m 

1000  m 

15. 

8.0 

19 

53.1202 

l/3(e) 

5  Mt 

0  m 

1500  m 

5.6 

3.5 

20 

53.1203 

1 /3(e) 

5  Mt 

0  m 

2000  m 

2.7 

2.0 

21 

53.6001 

1/4 

5  Mt 

0  m 

1000  m 

15. 

2.5 

22 

r  ^  r  r»r\o 

OO  •  Ou \JC 

1  M 

l/H 

5  Mt 

0  m 

1500  m 

4.1 

1.2 

23 

53.6003 

1/4 

5  Mt 

u  m 

2000  m 

1.8 

0.66 

24 

33.6011 

l/4(a) 

3  Mt 

U  m 

1000  m 

9.0 

1.4 

25 

33.6012 

l/4(a) 

3  Mt 

u  in 

1500  m 

2.6 

0.68 

26 

33.6013 

l/4(a) 

3  Mt 

2000  m 

1.1 

0.34 

27 

53.6011 

l/4(a) 

*  Mt 

1000  m 

15. 

2.2 

28 

53.6012 

l/4(a) 

5  Mt 

1500  m 

4.1 

1.1 

29 

53.6013 

1 /4(a) 

5  Mt 

2000  m 

1.8 

0.64 

30 

13.1211 

2/3 

1  Mt 

1000  m 

4.8 

2.2 

31 

13.1212 

2/3 

1  Mt 

0  m 

1500  m 

2.0 

1.1 

32 

13.1213 

2/3 

1  Mt 

2000  m 

1.0 

0.64 

33 

13.1012 

2/ 3(a) 

1  Mt 

0  m 

1500  m 

1.4 

0.84 
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APPROACH 


A  uniform  laydown  of  reentry  vehicles  (RV)  over  an  evenly  distributed  system 
of  targets  (shelters)  was  assumed.  The  ratio  of  RV  to  targets  was  used  to  des¬ 
cribe  the  attack  laydown  pattern.  Thus,  attack  scenario  1/2  refers  to  an  attack 
laydown  such  that  there  is  one  RV  attacking  each  pair  of  shelters.  The  suffixes 
on  the  various  scenarios  (e.g.,  1/2  (a))  refer  to  the  different  laydown  configura¬ 
tions  possible.  Using  this  nomenclature,  the  12  attack  scenarios  examined  were 
1 /2(a) ,  l/2(b) ,  l/2(c) ,  l/2(d),  1/3,  l/3(a),  l/3(c),  l/3(e),  1/4,  l/4(a),  2/3  and 
2/3(a)  (figures  are  in  Appendix  A).  Yields  of  1 ,  3  and  5  Mt  were  used.  Spacing 
between  shelters  was  1000,  1500,  or  2000  m.  Bursts  were  numbered  sequentially 
as  a  function  of  increasing  angle  clockwise  from  0°  (top  of  page)  and  were 
detonated  simultaneously  at  ground  level. 

Stations  were  not  colocated  with  shelters  (which  tend  to  fall  on  symmetry 
points)  for  the  following  reasons: 

1.  The  LAMB  code  underpredicts  peak  overpressure  at  symmetry  points. 

2.  Perfect  symmetry  will  not  occur  in  MX  construction. 

3.  Zero  CEP  is  highly  improbable. 

4.  Assuming  zero  CEP,  simultaneous  detonation  and  perfect  symmetry 
are  not  operationally  obtainable  environments. 

Thus,  measuring  overpressure-time  histories  at  nonsymmetry  points  provides  more 
realistic  waveforms  for  physically  plausable  attack  scenarios.  Stations  were 
placed  100  m  from  the  shelter  and  20°  off  a  direct  line  between  the  shelter  and 
the  nearest  burst  (Fig.  1).  Three  stations  were  required  to  define  overpressure 
and  overpressure  impulse-time  histories  for  the  1/3  and  l/3(a)  scenarios.  One 
station  sufficed  for  each  of  the  other  scenarios. 


station 

point  (Not  to  Scale) 


shelter  point 


nearest 

burst 


Figure  1.  Shelter-station  relation. 
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Test  cases  were  run  to  determine  the  number  of  bursts  to  be  used  in  each 
calculation.  These  test  cases  represent  the  worst  possible  attacks  (for  the 
scenarios  examined)  in  terms  of  the  number  of  bursts  reaching  the  station  point 
at  approximately  the  same  time.  The  results  of  the  test  cases  (Fig.  2  through 
7)  show  that  restrictions  on  the  number  of  bursts  in  a  calculation  could  be  made 
while  preserving  the  significant  overpressure-time  history  at  the  station  point. 
All  bursts  which  were  within  four  closest  burst-to-shelter  radii  of  the  shelter 
were  included  in  the  calculations  for  the  1/3  and  1/4  scenarios.  For  economy, 
the  number  of  bursts  was  restricted  to  12  for  each  of  the  other  scenarios. 

All  the  calculations  were  run  to  detonation  time,  plus  12  s.  All  of  the 
shock  waves  had  passed  over  each  station  by  this  time.  The  shock  front  velocities 
were  computed  in  the  following  manner: 

A  shock  velocity  function  subroutine,  VSHCK,  was  called  from 
the  LAMB  code  with  height  of  burst,  shock  radius,  height-of- 
station  point,  and  yield  as  input  parameters.  The  analytic 
expression  used  for  shock  velocity  was  obtained  from  the 
Rankine-Hugoniot  relation: 


U  = 


+ 


2y 


AP  W2 

P  / 


where 

U  =  shock  front  velocity 
c  =  ambient  speed  of  sound 
y  =  ratio  of  specific  heats  for  air 
AP  =  overpressure 
P  =  ambient  pressure 
Y  =  1.4  was  used  for  this  case 

The  direction  of  the  shock  was  determined  by  computing  the  angle  between  the 
burst  and  the  station  in  polar  coordinates. 
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RESULTS 

Overpressure  and  overpressure  impulse-time  histories  are  presented  in  Appendix 
B.  The  number  associated  with  each  peak  is  the  burst  number.  The  shock  arrival 
time  (s),  shock  velocity  (m/s),  shock  direction  (°),  burst  number,  and  burst  range 
(m)  are  presented  in  tabular  form.  The  attack  scenario  is  presented  at  each  over¬ 
pressure  plot  where  the  number  corresponds  to  a  burst  point  and  the  letter  s 
depicts  a  station  point  (Appendix  A).  Shock  direction  is  equivalent  to  the 
bearing  from  the  burst.  As  it  passes  over  the  station,  the  shock  is  travelling 
parallel  to  the  bearing  from  burst  station.  Burst  range  refers  to  the  distance 
between  the  burst  and  station. 
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Figure  2.  Overpressure  versus  time-2/3(a)  attack  (with  added  bursts) 
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Figure  3,  Overpressure  pulse  versus  time-2/3(a)  attack  (with  added  bursts). 
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Figure  4.  Overpressure  versus  time-2/3  attack  (with  added  bursts). 
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Figure  5.  Overpressure  impulse  versus  time-2/3  attack  (with  added  bursts) 
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Figure  6.  Overpressure  versus  time-l/2(c)  attack  (with  added  bursts). 
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Figure  7.  Overpressure  impulse  versus  time-l/2(c)  attack  (with  added  bursts) 
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APPENDIX  A 

The  following  12  figures  illustrate  the  various  attack  scenarios  used  in  this 
multiburst  study.  Included  in  these  figures  are  launch  points,  aim  points,  calcu- 
lational  burst  points  and  station  points.  The  various  attack  scenarios  are 
described  as  a  ratio  of  RV  to  targets.  For  example,  a  1/2  scenario  refers  to 
one  RV  attacking  each  pair  of  targets. 
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APPENDIX  B 

The  following  figures  are  overpressure  and  overpressure  impulse-time  histories 
of  the  various  attacks  described  in  Appendix  A.  The  burst  number  associated  with 
the  overpressure  peaks  corresponds  to  the  attack  scenario  burst  number  inset  on 
each  plot.  In  addition,  the  shock  arrival  time,  shock  velocity,  shock  direction, 
burst  number  and  burst  range  are  presented  on  the  overpressure  plots. 
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